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ABSTRACT 


A  cccputer  model  was  developed  to  simulate  a 
projectile  ram- propellant  feed  system  for  a  Liguid 
Propellant  Gun.  Using  a  lumped  parameter  approach,  a 
set  of  sinultaneous  differential  eguations  was  derived 
for  tbe  complex  interaction  of  the  propellant  fluid, 
the  driving  injector  and  the  projectile.  The  computer 
model  was  verified  against  a  20  mm  experimental 
apparatus.  Injector  displacement,  projectile 
displacement,  and  chamber  pressure  were  compared  for  a 
nominal  driving  pressure  of  140  psi.  The  important 
systen  parameters  affecting  projectile  ram  time  and 
chamber  pressure  oscillations  were  investigated  and 
potential  prctlem  areas  for  testing  with  actual 
propellant  were  identified. 
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INTRODUCTION 


As  the  state  cf  the  art  of  weapons  technology  improves, 
mere  complex  systems  are  continually  being  developed.  To 
meet  future  needs  these  systems  must  have  improved 
performance,  less  weight  and  volume,  less  complicated 
logistics,  and  at  the  same  time  be  more  cost  effective  than 
current  sjstems.  One  candidate  for  these  future  weapon 
systems  is  the  liquid  propellant  gun  (LPG) .  Numerous 
investigaticcs  over  the  past  thirty  years  have  demonstrated 
its  potential  if  the  technology  can  be  developed  to  utilize 
a  liguid  propellant  in  a  large  caliber  gun. 

One  important  area  that  needs  investigation  is  the 
aspect  of  fluid  prcpellant  handling.  The  ordnance  designer 
must  knew  the  necessary  propellant  flow  rates  and  pressure 
regimes  needed  to  achieve  the  desired  system  performance. 
These  parameters  must  be  within  allowable  safety  margins  for 
the  handling  and  use  of  the  explosive  propellant.  One 
possible  design  of  an  LPG  which  would  allow  a  very  high 
firing  rate  utilizes  the  propellant  to  ram  the  loaded 
projectile.  Combining  the  normally  separate  propellant  load 
and  ram  cycles  decreases  the  reguired  novement  of  mechanical 
parts.  Thus,  using  this  method,  it  is  feasible  that  firing 
rates  in  excess  of  fpur  to  five  times  present  rates  could  be 
achieved.  This  thesis  was  directed  toward  identifying  the 
important  fluid  dynamic  parameters  involved  with  such  a 
prcjectile  ram  feed  system. 


II.   BACKGBOUNC 


During  1S74  and  early  1975,  an  investigation  was 
conducted  at  the  Naval  Postgraduate  School  to  study,  both 
analytically  and  experimentally,  the  fluid  dynamics  of  a 
liguid  propellant  under  conditions  similar  to  those  which 
would  exist  in  a  rapid-fire  LPG  feed  system.  The  results  of 
this  investigation  were  to  b€  used  to  establish  such  LPG 
design  and  performance  parameters  as  time-to-load,  injection 
supply  pressure,  injection  system  configuration,  ullage, 
charge-tc-mass  ratio,  caliber  size,  and  projectile  mass.  It 
was  hoped  that  this  investigation  would  identify  any 
potential  problem  areas  for  further  detailed  research. 

The  basic  objective  of  the  experimental  portion  cf  the 
investigation  was  to  identify  what  fluid  dynamic 
characteristics  of  a  liguid  propellant  feed  system  would 
limit  leading  times  and  hence  rates  of  fire.  A  20  mm 
experimental  model  of  a  basic  propellant  feed  system  was 
designed  and  built.  Data  on  injector  displacement,  breech 
chamber  pressure,  and  ram  gas  pressure  (input  driving 
pressure) ,  were  recorded  for  driving  pressures  between  50 
and  220  psig  in  10  psi  increments.  As  reported  in  Ref.  1, 
it  was  found  that  the  instantaneous  behavior  of  the  chamber 
pressure  was  the  result  of  a  complex  interaction  of  inertia 
fcrces,  visccus  forces,  and  the  unsteady  motion  of  the 
fluid.  Ike  experiments  demonstrated  that  frictional  and 
inertial  effects  were  significant  during  the  movements  of 
the  injector  and  the  projectile  slug.  Once  the  projectile 
slug  stopped,  the  effect  of  entrapped  gas  in  the  fluid 
caused  large  breech  chamber  pressure  oscillations.  On 
several  runs,   sub-atmospheric   pressures   were  experienced 


which  suggested  the  possibility  of  cavitation  and  hence 
vapor-phase  igniticn-  Over  all,  it  was  found  that  the  ram 
tine  had  a  guadratic  dependence  upon  ram  pressure. 

The  analytical  portion  of  the  investigation  was  directed 
toward  predicting  the  pressure  and  flew  rate  of  the 
propellant  atd  the  projectile  slug  motion  during  an  LPG 
loading  cycle.  The  mathematical  model  which  was  developed 
and  programmed  on  an  analog  computer  predicted  the  position, 
velocity,  and  acceleration  of  the  projectile  anc  the 
pressure  at  various  points  in  the  system  as  functions  of 
time.  The  mcdel  was  tested  against  experimental  results  and 
found  to  he  adequate  for  the  prediciton  of  projectile  ram 
time.  An  analysis  was  performed  using  the  model  to  indicate 
the  areas  cf  system  redesign  likely  to  be  most  profitable 
and  to  cttain  preliminary  predictions  of  LEG  loading  system 
performance  under  a  variety  of  design  conditions.  The 
results  cf  the  analytical  study,  as  well  as  a  summary  cf  the 
experimental  work  is  contained  in  Bef.  2. 

The  aforementioned  model  was  the  starting  point  fcr  the 
present  study.  To  better  understand  the  derivations  in  the 
fcllowing  sectiens,  Figure  (1)  indicates  the  geometry  cf  a 
hasic  prc:ectile  ram  propellant  feed  system.  An  injector 
chamber  is  filled  with  propellant,  which  is  then  pumped  into 
the  gun  treech  by  applying  ram  pressure  to  the  ram  side  of 
the  injector  piston.  This  can  be  accomplished  by  using  high 
pressure  gas  from  an  accumulator  or  by  using  an  hydraulic 
drive  systen.  The  force  exerted  by  the  injector  piston 
drives  the  fluid  propellant  through  the  connecting  line  into 
the  gun  treech.  The  rapidly  accelerated  propellant  drives 
the  already  leaded  projectile  to  its  seated  position  ready 
for  firing. 
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Figure    1  -   EASIC  GEOMETBY 
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The  ccuplex  nature  of  this  dynamic  ram  process  is  not 
easy  to  describe  with  mathematical  equations.  Many 
researchers  have  approached  the  problem  using  wave  mechanics 
and  partial  differential  equations.  These  equations  were 
then  modified  according  to  the  propagation  characteristics 
of  the  system  and  the  various  boundary  conditions 
encountered  therein.  This  usually  led  to  involved  finite 
difference  methods  of  solution  on  a  digital  computer. 

Another  possihle  approach  considers  the  kinetic  energy 
involved  in  tbe  ran  process.  This  is  the  approach  that  was 
used  to  derive  the  governing  eguations  for  the  previously 
described  analog  computer  model.  Beginning  with  the  input 
side,  a  force  balance  on  the  injector  piston  yielded  an 
eguaticn  fcr  the  pressure  in  the  injector  chamber  in  terms 
of  the  injector  motion.  Writing  Bernoulli's  eguations  for 
head  losses  at  the  expansion  to  the  breech  chamber,  the 
ccntracticn  tc  the  connecting  line,  and  an  orifice  in  the 
connecting  line  resulted  in  an  equation  fcr  the  pressure 
drcp  in  the  connecting  line  in  terms  of  the  square  cf  the 
fluid  velocity  and  the  fluid  acceleration.  A  force  balance 
on  the  projectile  slug  yielded  a  third  equation  which 
described  the  breech  chamber  pressure  in  terms  of  the  slug 
motion.  Tc  solve  these  equations,  seme  method  of  relating 
the  motions  in  the  injector,  connecting  line,  and  treech 
chamber  had  to  te  established.  It  was  assumed  that  the 
fluid  was  incompressible  and  therefore  these  motions  were 
identical.  To  facilitate  the  analog  simulation  it  was  also 
assumed  that  the  irput  ram  pressure  was  a  step  input.  In 
the  process  of  converting  these  eguations  to  an  integral 
fcrm  for  wiring  on  an  analog  computer,  all  variables  were 
normalized  tc  permit  scaling. 

Despite  these  limiting  assumptions,  the  analog  computer 
model  was  able  to  predict  the  time  to  ram  the  projectile  and 
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the  chanter  jressures  which  cccur  during  the  time  that  the 
injector  is  in  motion.  Due  to  the  incompressitility 
assumption,  the  model  was  unable  to  account  fcr  the  pressure 
transients  cr  "water  hammer"  type  pressure  oscillations 
which  occurred  when  the  projectile  was  seated  and  the 
propellant  fluid  was  being  decelerated.  These  pressure 
oscillations  are  important  to  the  system  designer  because 
the  peak  breech  chamber  load  pressure  is  experienced  during 
these  oscillations.  It  is  also  possible  that  these 
oscillations  could  interfere  with  the  uniformity  of  the 
ignition  and  subseguent  combustion  of  the  propellant. 

A  straight  fcr.ward  attempt  was  made  to  extend  this 
analog  model  ry  first  converting  it  to  a  digital  computer 
program.  This  allpwed  an  accurate  representation  of  the 
input  ram  pressure  and  the  use  of  unnormalized  variables  to 
be  incorporated.  These  improvements  increased  the  accuracy 
of  the  model  rut  the  effect  of  the  compressibility  of  the 
propellant  fluid  still  was  not  taken  into  account.  Trying 
to  add  compressibility  effects  to  this  model  by  adding  time 
derivatives  cf  the  pressure  terms  to  the  governing  eguations 
unlinked  the  motions  of  the  injector,  the  connecting  line 
fluid,  and  the  projectile  slug  causing  a  problem  with  too 
many  unknown  variables  for  the  number  of  eguations  involved. 

It  was  decided  that  a  new  approach  should  be  tried  which 
would  adeguately  describe  the  system  pressure  oscillations. 
The  approach  would  feature  an  engineering  model  which  would 
be  easy  to  use,  adaptable  to  any  LPG  feed  system,  and  not 
obscure  the  interaction  of  system  parameters  by  complicated 
mathematics. 
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III.   DERIVATION  CF  COMPUTER  MODEL 


LUMPEI  I2EAMETEB  APPROACH 


The  approach  that  showed  the  greatest  promise  for 
modeling  the  LPG  feed  system  was  that  of  the  fluid 
transmission  line  cpncept.  This  approach,  which  has  recome 
popular  in  the  last  five  to  ten  years,  is  based  on  a 
pressure-vcltage  and  flow  velocity-current  analogy  with 
Electrical  Engineering  determinations.  It  is  an  outgrowth 
of  the  large  amount  of  effort  that  has  been  devoted  to 
investigating  fluid  line  transients.  Reference  3  is  a  good 
survey  of  this  field. 

For  complex  systems  this  approach  is  usually  simplified 
fcy  using  the  approximation  of  lumped  parameters.  The 
effects  of  fluid  inertia,  capacitance,  and  resistance  are 
"lumped"  and  considered  to  act  only  in  discreet  areas  cf  the 
system.  This  results  in  a  reduction  cf  the  unknown 
parameters  due  tc  the  assumed  lack  of  interaction  cf  the 
different  fluid  effects.  This  approximation,  however, 
results  in  the  .necessity  of  using  several  empirical 
constants  which  must  be  determined  by  fitting  experinental 
data.  The  ordinary  differential  eguations  that  are  derived 
frcm  this  method  can  be  solved  either  by  Laplace 
transf crmaticn  or  fcy  computer  integration. . 

To  consider  the  effect  of  fluid  inertia  it  is  assumed 
that  cnly  pressure  and  inertia  forces  are  present  and  that 
compressibility  effects  are  negligible  in  the   volume   under 
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consideraticr.   Then, 

ay 

pl-p« -/>«■<? 

where  C        is   the   propellant  fluid   density,   L    is   a 

characteristic   length  and   V  =V  =V   because   the   flew   is 

1   2 

considered  incompressible  for  this  building  block  (Fig.2A). 


To   consider   the  effect  of   fluid   capacitance,  it  is 

assumed  that  only  compressibility  effects  are  important,  and 

that   inertia  and  resistance  effects  may  be  neglected  in  the 

vcluire   UDder   consideration.      Therefore, 

L  dP 

V  -  V  = 

1     2   K  dT 
where   K   is   the   effective   bulk  modulus  of  the  propellant 

fluid,  I  is  a  characteristic  length  and  P  =F  =P  (Fig.  2B) . 

1   2 


Eecause  fluid  resistance  can  be  affected  by  sc  many 
different  parametexs,  it  is  impossible  to  write  a  general 
eguation  describing  the  pressure  drop  due  to  fluid 
resistance.  It  is  best  to  treat  it  empirically  using  an 
experimentally  derived  figure.  Thus,  in  the  volume  under 
consideraticr, 

P  -  P  =  Rv  V 

1     2       . 

where   Ev   is  a  function  of  fluid  velocity  and  V  =V  =V  (Fig. 
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2C) .   Changes  in  cross  sectional  area  can  be  accounted  for 
by  using  appropriate  area  ratios. 

These  results  specify  three  building  blccks  which  can  be 
ccibined  in  any  sequence  to  model  the  dynamic 
characteristics  of  a  system.  The  complexity  of  the  model 
can  be  increased  to  any  degree  necessary  by  including  more 
and  more  cembinatiens  of  these  three  basic  building   blocks. 
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A)    Inertia 


P    -  P 


dv 


B)     Capacitance 


V  =  V  =  V 
i   2 


L  dp 
y  -  y  .  ._ 

1     2   K  dT 


p  =  p  =  p 

1     2 


C)  Resistance 


P  -  P  =  Rv  V 

1     2 


V  =  V  =  V 

1     2 


Eigur€  2  -   LOMfED  EARAMETEB  MCEULES 
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I-PG  SYSTEM  GOVERNING  EQUATIONS 


Several  different  combinations  of  building  blocks  were 
tried  iD  attempting  to  model  the  experimental  LPG  system. 
The  experimental  results  showed  that  inertial  and  resistance 
blocks  should  te  included  in  the  model.  The  fact  that 
pressure  transients  pccurred  in  the  breech  chamber  indicated 
that  a  capacitance  block  should  also  be  included.  To 
ninimize  the  use  of  computer  time,  it  was  decided  to  start 
with  the  mcst  simple  model  thought  adequate  and  work  toward 
more  complex  models  as  the  guality  of  the  computer  results 
dictated. 

The  first  attempt  considered  that  ccmpressibilty  effects 
would  dcttinate  in  the  injector  chamber  as  the  accelerating 
piston  interacted  liith  the  propellant  fluid.  The  resistance 
effects  cf  changes  in  the  feedline  cross  sectional  diameter 
and  associated  orifices  were  modeled  as  one  resistance  block 
at  the  beginning  of  the  connecting  line.  Inertial  effects 
*ere  thought  to  dominate  in  the  connecting  line.  This  would 
be  particularly  true  in  large  scale  models  where  the  length 
of  the  ccnnecting  line  would  be  very  large.  It  was  thought 
that  compressibility  effects  should  again  dominate  in  the 
breech  chamber  as  the  accelerating  propellant  fluid  drove 
the  projectile  dcwn  the  breech  chamber.  Results  in  this 
case  were  governed  by  inertial  effects  because  the 
calculated  ccnnecting  line  velocity  was  too  great.  This  led 
to  a  decrease  in  the  chamber  pressure  as  the  projectile 
accelerated  until  a  large  negative  chamber  pressure  existed. 
In  addition  to  this  deficiency,  the  computed  system  pressure 
cscillaticns  were  not  sufficiently  damped. 

To  try  tc  adjust  the  connecting  line   velocity,   another 
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resistance  tlcck  was  added.  As  will  be  seen  in  the  section 
discussing  the  parameters  which  affect  the  pressure 
oscillations,  the  \alue  of  this  resistance  coefficient  is 
proportional  to  the  amount  of  damping  in  the  pressure 
oscillations .  Even  with  a  value  which  gave  the  proper 
danping,  the  breech  chamber  pressure  was  dominated  ry  the 
pressure  loss  due  to  the  accelerating  fluid  in  the 
connecting  line. 

To  reduce  the  dominance  of  the  fluid  inertia  effects,  it 
was  decided  tc  restructure  the  model.  To  model  the  injector 
chamber  it  *as  decided  to  consider  the  effect  of  fluid 
capacitance  and  inertia.  The  connecting  line  was  modeled  by 
t*c  resistance  blocks  separated  by  a  fluid  capacitance 
block.  The  breech  chamber  was  modeled  using  a  fluid 
capacitance  and  a  fluid  inertia  block,  as  was  done  for  the 
injector  chamber.  This  mcdel  gave  the  final  results 
described  in  the  section  comparing  computer  and  experinental 
data.  Ey  combining  these  building  blocks,  the  governing 
eguatiens  became:   (See  Fig.  (1)  for  notation) 

For  the  injectcr: 

1.    Pr  A    -  P  A    =  Mp  V    +  Kfs  V    +  Pdi  A 
ill  1  li 


2.    V    -V    -  J-  P, 

1         2       K 


P    - P    =  pl_    V 
2  3      r    1      2 


For    the    connecting   line: 
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u.  P  -  P  =  Rv  V 

3    4     14 


L2  • 
5.  V   -  V  =  —  P 
4    5    K   4 


6.  p   -  p  =  RV  V 
5    6     2   5 


For  the  rreech  chamber: 


7.   P   -  P   =  pi      V 
6    7    ~  3   6 


L3  ' 
8.  V   -  V  =  -  P 
7    8   K   8 


• 


P  A  -  Pds  A  =  Ms  V  +  Kfs  V 

8   3,3       8        8 


These  equations  were  subsequently  modified  tc  account 
for  changes  in  cross-sectional  area  and  changing  geometry  as 
the  injector  chamber  decreased  in  volume  and  the  fcreech 
chamber  increased.  Once  the  projectile  reached  the  end  of 
the  breech  chamber,  the  force  balance  equation  was  no  longer 
considered  applicable.  At  this  time,  it  was  noted  that  the 
renaining  equations  for  the  connecting  line  and  the  fcreech 
chamber  could  be  combined  to  form  a  second  order 
differential  equation  for  the  breech  chamber  pressure  which 
described  tie  "water  hammer"  type  pressure  oscillations 
displayed  by  the  experimental  data.   Thus, 
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«   A3  Rv2  . 
10.  P   +  r2  -f-  P 


A2  />L3   8    /)L3L3   8   ^31.3   4 

Onc€  the  injector  piston  reached  the  end  of  its  travel, 
the  injectcr  chamber  no  longer  existed  acd  the  governing 
eguaticns  were  considered  as  no  lcnger  applicable. 
Ultimately,  the  injector  chamber  pressure  was  propagated 
doiin  the  connecting  line  as  the  propellant  fluid  came  to 
rest.  Throughout  these  analyses,  input  ram  pressure  was 
modeled  using  expcrential  terms  to  fit  experimental  data. 

The  eguaticns  from  the  lumped  parameter  approach  fcr  the 
prcpellant  fluid  were  combined  with  the  force  balance 
eguations  fcr  the  injector  piston  and  the  projectile  using 
State  Variable  methods.  Subsequently,  a  computer  program 
was  constructed  to  solve  these  state  variable  relations 
using  a  fotrth  crder  fiunga-Kutta  integration  routine  for 
sinultanecus  first  order  differential  eguations  which  was 
developed  at  the  Naval  Postgraduate  School.  A  brief 
description  of  the  state  variable  method  and  a  listing  of 
the  computer  program  can  be  found  in  Appendix  A. 

The  following  is  a  listing  of  the  nomenclature  and  the 
values  of  constants  used  in  the  derivation  of  the  governing 
eguations  anc  the  computer  program: 

A1    Injectcr 

Cross-Sectional  Area 

A2   Connecting  Line 

Cross-Sectional  Area 

A3    Ereech  Chamber 

Cross-Sectional  Area 

IB   Earn  listen 

Cross-Sectional  Area 
11    Injector  Length 
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1.77 

2 
in 

11.42 

2 

Cff 

.255 

2 
in 

1.65 

2 

Cff 

.49 

2 
in 

3.16 

2 

Cff 

1.77 

2 
in 

11.42 

2 

Cff 

1.60 

in 

4.22 

Cff 

12 
L3 

ME 

MS 

KFP 

KFS 

FDI 
ECS 

K 

EHC 

PB 
P1 
V1 
V2 
F4 
V6 
E8 
V8 


Connecting  Lire  Length    16.0 
Breech  Chamber  Length    5.0 


.0052 


Injectcr  Piston  Mass 

Projectile  Mass 

Injectcr  Frictipn 

Factcr 

Projectile  Friction 

Factcr 

Injector  Eack  Pressure 

Projectile  Back 

Pressure 

Effective  Bulk  Modulus 

Prcpellant  Density 


Bam  Pressure 

Injectcr  Chamber  Pressure 
Injectcr  Pistcn  Velocity 
Injectcr  Chamber  Exit  Velocity 
Connecting  Line  Pressure 
Breech  Chamber  Entrance  Velocity 
Ereech  Chamber  pressure 
Projectile  Velocity 


in       40.6   ci 

in       12.7   ce 

2 
lbrsec.   913.8  gn 

in 


.00053 

lb-sec 
in 

93.1 

gn 

.001 

lb-sec 
in 

175.7 

gn/sec 

.001 

lb-sec 
in 

175.7 

ga/sec 

24.0 

PSI 

163.3 

kEa 

.50 

PSI 

3.4 

kEa 

3200 

PSI 

6808 

kEa 

.000093 

2 

lb-sec 

.994 

3 

gn/cm 

in 


Computed 
Computed 
Computed 
Computed 
Computed 
Computed 
Computed 
Computed 


The  treatnent  of  several  areas  of  the  computer  model 
continually  reappeared  as  reguiring  refinement.  The 
governing  differential  equations  worked  uell  during  the 
dynamic  portion  cf  the  feed  cycle  but  experienced 
difficulties  during  the  initial  and  final  static  periods. 
The  handling  cf  static  friction  and  back  pressure  as  a 
ccrstant  value  created  the  possibility  of  negative 
velocities  urtil  the  driving  ram  pressure  overcame  the 
systen   tack  pressure.   These  negative  values  never  occurred 
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in  the  real  system  due  to  the  gecmetric  restraints  on  the 
injector  piston  and  projectile;  therefore,  the  coaputer 
model  had  to  he  manipulated  to  maintain  this  condition. 
Unfortunately  the  need  for  simultaneous  solution  cf  the 
governing  eguations  made  it  difficult  to  manipulate  the 
initial  conditions  without  greatly  increasing  the  complexity 
of  the  computer  model.  After  the  projectile  stopped,  the 
transiticr  frcm  the  force  balance  equation  en  the  projectile 
to  an  eguaticn  describing  the  fluid  velocity  in  response  to 
the  pressure  transients  was  awkward.  No  simple  differential 
eguation  describes  the  complexity  of  the  wave  mechanics 
involved  with  the  reflection  of  the  pressure  waves  in  the 
system.  In  this  regard  several  alternate  approaches 
involving  seguential  alterations  of  the  governing  eguations 
were  tried  with  varying  degrees  of  success. 
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I*.   DESCBIPTION  CF  EXPERIMENTAL  APPARATUS 


In  crdei  that  the  comparison  of  the  experimental  and 
ccuputer  generated  data  can  be  fully  understood,  a  brief 
description  cf  the  NPS  experimental  apparatus  and  the 
conduct  cf  tie  associated  experiments  is  included  (See  Ref. 
1,  pages  17-27). 

The  test  chamber  was  fabricated  frcm  a  three-inch  CD. 
Lucite  cjlirder,  18  inches  long,  bored  to  a  20  nm  inside 
diameter  and  fitted  with  aluninuu  end  caps.  The  chamber  was 
loaded  with  a  trass  slug  weighing  93  grams  which  rode  en  two 
graphite  filled  Teflpn  sealing  rings.  The  brass  slug,  which 
sinulated  the  projectile,  was  cycled  from  the  breech  end  of 
the  chamber  tc  the  barrel  end  and  returned  to  the  treech 
end,  completing  one  hypothetical  firing  cycle. 

Because  cf  the  desire  to  vary  the  charge  to  mass  ratio, 
a  variable  chamber  velocity  was  necessary.  To  accomplish 
this  with  one  chamber,  a  volume  control  retaining  rod  was 
designed  intc  the  system.  This  brass  rod,  bored  to  allow 
gas  to  pass  its  length,  was  threaded  through  a  plate  which 
was  attached  to  the  barrel  end  cap  holding  the  red  in  the 
chamber.  The  rod,  wiiich  has  a  Teflon  disc  en  the  end,  not 
enly  established  the  volume  of  the  test  chamber,  but 
provided  a  buffer  stop  for  the  slug  at  the  end  of  its 
forward  mcticn.  Another  Teflon  buffer  was  affixed  to  the 
breech  eod  cap  to  cushion  the  slug  in  return  motion. 

The  LEG  simulator  started  a  simulated  firing  cycle  with 
the  slug  at  the  breech  end  of  the  empty  chamber,  as  shewn  in 
Fig.  (3). 
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This  is  the  ready-to-ram  position.  The  simulated 
prcpellant,  distilled  water,  was  then  introduced,  ramming 
the  slug  tc  the  opposite  (barrel)  end  of  the  chamber  as  the 
chamber  was  filled.  This  was  accomplished  by  applying  gas 
pressure  tc  a  power  piston  which  drove  the  injector  piston. 
The  injector  piston  forced  the  propellant  past  a  flow  check 
valve  and  into  the  chamber.  This  placed  the  slug  in  the 
ready-tc-fire  cr  ic-battery  position.  In  an  actual  gun,  the 
prcpellant  would  re  ignited  at  this  time  in  the  cycle.  Due 
to  laboratory  constraints,  an  expulsion  system  was  used. 

An  IVDI  displacement  transducer  was  manufactured  and 
mounted  next  tc  the  injector  piston.  The  LVDT  was  attached 
tc  the  connecting  rod,  between  the  power  piston  and  the 
injectcr  pistcn.  The  volume  of  the  liguid  being  placed 
under  pressure  during  each  shot  was  measured  by  filling  the 
system  in  the  ready-to-ram  position  and  then  draining  it 
intc  a  graduated  beaker.  By  measuring  the  displacement  of 
the  injectcr  piston  head  with  the  LVDT,  the  volumetric  rate 
of  fluid  injection  during  the  ram  stroke  was  obtained. 

Two  pressure  taps  were  drilled  in  the  test  chamber. 
These  taps  were  located  as  close  to  the  breech  end  as 
possible,  cne  at  2CT  degrees  frcm  top  center  and  the  other  at 
20  degrees  from  bottom  center  in  a  counterclockwise 
direction,  as  seen  from  the  breech  end  of  the  chamber.  Only 
the  bottom  location  was  used  to  record  data. 

A  4-chanrel  Hewlett-Packard  3960  Magnetic  Tape  Recorder 
was  used  to  1M  record  desired  data  during  system  operation. 
For  each  cata  run  signals  were  recorded  from  two  Kaman 
diaphragms  type  ,(1000  psi)  pressure  transducers,  one 
connected  to  the  breech  pressure  tap  on  the  test  chamber, 
and  the  ether  to  the  gas  injection  side  of  the  power  piston. 
The  pressure  signals  were  processed  with  a  Kaman  Digi-Vit 
Eeadout  Dnit  which  also  provided  a  visual  (digital)  display. 
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Data  fron  the  IVDT  displacement  transducer  was  also  recorded 
en  the  tape. 

A  Brush  Recorder  (Mark  280)  was  used  to  obtain  a  visual 
display  cf  the  recorded  data.  By  transcribing  the  desired 
signals  en  the  Magnetic  Tape  Becorder  at  a  tape  speed  cf  15 
feet  per  minute  and  playing  them  back  into  the  Brush 
Becorder  at  5-3/4  feet  per  minute,  the  time  scale  of  the 
output  was  expanded  by  a  factor  of  four  on  the  Brush 
recordings  (viz.,  from  a  real-time  maximum  cf  200  mm/sec  to 
a  delayed  tiie  maximum  of  800  mm/sec.) 
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CCKPjURISON  op  COMPOTEB  AND  EXPERIMENTAL  DA1A 


As  menticned  previously,  experimental  data  was  taken  for 
ran  gas  pressures  ranging  from  50  to  220  psig.  Fron  this 
collection  of  data,  one  run  at  140  psig  was  selected  as 
representative  of  system  performance.  It  was  felt  that  any 
model  which  *culd  adeguately  describe  systeir  operation  at 
this  intermediate  pressure  would  be  valid  for  the  entire 
range  of  expected  1PG  driving  pressures. 

Figure  (U)  is  a  comparison  of  computed  and  recorded  data 
for  injector  pistot  displacement.  As  can  be  seen,  the 
agreement  is  very  good.  This  is  not  very  significant  in 
that  all  iicdels  tried,  as  well  as  the  original  analog 
computer  mocel,  were  able  to  correctly  predict  injector 
piston  metier. 

Figure  (5)  is  a  comparison  of  analytical  results  and 
experimental  data  for  breech  chamber  pressure.  The  computer 
model  fellows  the  shape  of  the  experimental  curve  fcr  the 
duration  of  the  tine  that  the  projectile  is  in  motion  (0-30 
msec) .  It  oscillates  very  rapidly  but  does  not  fall  off  to 
a  2erc  value  at  the  end  of  the  projectile  moticn.  The 
freguency  of  these  initial  oscillations  is  approximately  1.0 
kh2  which  is  within  the  -frequency  range  of  the  Kaman 
pressure  transducers  used  in  testing.  However,  the  mounting 
of  these  transducers  within  a  connecting  cavity  instead  of 
flush  with  the  chamber  wall  could  have  led  to  them 
experiencing  a  reduced,  lagging  frequency  response.  It  is 
felt  that  seme,  if  not  all,  of  the  computed  pressure 
oscillations  must  exist  as  evidenced  by  the  close 
correspondence   of   the   first  two   peaks  in  Fig.  (5) .   The 
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accuracy  cf  the  experimental  pressure  reading  could  also  be 
considered  as  being  reduced  by  the  location  of  the  pressure 
transducer  ic  the  test  setup.  It  is  possible  that  the 
monitcricg  pressure  tap  was  too  close  to  the  end  cf  the 
breech  chamber  and  may  not  have  sensed  the  full  chamber 
pressure  during  the  latter  part  of  the  projectile^  travel 
when  the  fluid  velocity  is  greatest.  It  is  felt  that  to 
adeguately  cescrire  the  pressure  decrease  toward  zero,  it 
would  be  necessary  to  completely  account  for  exact  changes 
in  system  cecmetry.  To  do  this  would  require  using  a 
distributed  parameter  approach  featuring  variable  lengths 
fcr  the  in^ectcr  chamber  and  breech  chamber. 

Ihe  presscre  transients  which  occur  after  the  projectile 
steps  are  ccnplex  interference  phenomena  which  are  not  fully 
described  by  the  mcdel.  However,  the  peak  pressure  and  the 
natural  freguency  of  oscillation  of  the  computer  model 
cempare  fevcrably  with  the  experimental  data. 
Onfortunatelj,  the  damping  characteristics  cf  the  model  do 
not  fcllcM  the  experimental  data  well.  The  system 
parameters  which  determine  these  values  will  be  discussed  in 
the  next  section. 

It  shculd  be  reiterated  that  the  lumped  parameter 
approach  is  ar  approximation,.  In  a  distributed  parameter 
approach,  pressure  and  velocity  would  vary  ccntinuouslj  with 
distance  within  the  system  as  well  as  with  time.  In  the 
lumped  paraneter  approach,  these  changes  are  assumed  to 
occur  onlj  at  the  input  and  output  of  a  building  block; 
therefore,  pressure  and  velccity  are  considered  constant 
within  the  building  block.  This  will  always  lead  to  some 
discrepancies  when  comparing  model  data  to  experimental  data 
taken  at  at  a  fixed  point. 

Figure  (6)  shows  the  simultaneous  pressure  history  for 
several  areas   of   the  LPG   feed  system.    The   input  ram 
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pressure,  injector  chamber  pressure,  connecting  line 
pressure,  ard  breech  chamber  pressure  are  graphed  as 
functions  cf  time. 

Figures  (7)  and  (8)  show  the  simultaneous  velocities  at 
several  pcirts  as  predicted  by  the  computer  model.  The 
irjectcr  pistcE  velocity,  fluid  propellant  velocities  at  the 
beginning  ard  end  of  the  connecting  liDe,  and  projectile 
velocity  are  plotted  as  functions  of  time. 

The  injectcr  and  projectile  displacements  as  functions 
of  time  are  included  as  Figure  (9) .  As  can  be  seen  from 
Figures  (7) ,  (8) r  and  (9) ,  the  movement  of  the  injector 
piston  is  much  mere  rapid  than  that  of  the  projectile.  It 
is  felt  that  this  velocity  difference  is  very  nuch  a 
fuiction  cf  system  geometry  and  should  not  be  considered  as 
a  generali2ation  for  all  LPG  feed  systems.  The  NPS 
experimental  apparatus  being  modeled  had  only  small  pressure 
drops  between  injector  and  breech  chambers.  In  addition, 
the  projectile  slug  was  an  order  of  magnitude  lighter  than 
the  injectcr  piston.  In  large  scale  systems,  with 
significant  pressure  drops  and  very  massive  projectiles,  it 
is  guite  possible  that  the  injector  piston's  full  stroke 
will  occur  significantly  before  the  projectile  has  seated. 
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VI.   ANALYSIS  OP  SYSTEM  PERFORMANCE 


PARABETEIS  AFFECTING  PROJECTILE  RAM  TIME 


Several  parameters  can  be  varied  to  "tune"  the  model  to 
correct  ram  time.  In  equations  (1)  and  (9)  (the  force 
balances  en  the  injector  piston  and  projectile)  empirical 
constants  are  included  which  account  for  any  static  friction 
and  back  pressure  (PDI  =  injector  back  pressure,  PDS  = 
prcjectile  rack  pressure).  In  equations  (4)  and  (6), 
resistance  coefficients  are  introduced  to  account  for  system 
resistances.  Finally  in  equations  (2),  (5),  and  (8), 
prcpellant  effective  bulk  modulus  values  appear. 

The  value  for  the  injector  and  projectile  total  back 
pressure  is  difficult  to  determine  precisely.  Reference  1 
cited  an  experimentally  determined  value  of  24  psi.  The 
analog  computer  model  presented  in  Ref.  2  used  a  value  of  46 
psi  for  a  ran  pressure  of  140  psi.  The  present  model  uses  a 
value  of  24  psi  for  PDI  and  0.5  psi  for  PDS,  or  a  total  of 
24.5  psi  tack  pressure,  which  agrees  with  the  experimental 
value.  ihe  value  of  system  back  pressure  will  vary  with 
driving  pressure  and  the  configuration  of  each  LPG  system's 
injector  and  breech  chambers. 

Some  tasic  research  has  been  conducted  to  relate  the 
value  of  tie  resistance  coefficients  to  fluid  properties. 
Unf ortunatelj  these  studies  which  are  described  in  Refs.  4 
and  5,  dealt  with  steady  state  fluid  flow  in  cocstant 
diameter  lines   without   flew   restrictions,   making   their 
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results  inapplicable  to  the  present  LPG  feed  system  prcblem. 

The  value  cf  EV  can  be  related  to  the  damping  ratio  of   the 
2 

"water   hammer"  pressure  oscillations.   The  value  of  RV   was 

determined  tj  an  iterative  process  to  obtain  the  best  fit  to 

experimental  data.  Figure  (10)  shows  the  dependence  of 
ram-time   en   BV   and  RV  ,   keeping  all   ether   variables 

constant.    Is   expected  intuitively,   increasing  the  fluid 

resistance  creates  a  larger  system  pressure  drop,  resulting 
in  a  lewer  pressure  exerted  on  the  face  cf  the  projectile 
and  hence  slcwer  ram  times.  The  slope  cf  the  two  curves  are 
almost  identical  sp  that  no  significant  advantage  would  be 
achieved  ty  system  designs  which  try  to  minimize  either  RV 

or   RV   to  the  detriment  of  the  other.   It  should  be  noted, 
2 

however,  that  higher  values  for  RV   do  tend  to  slow  down  the 

ran  time  mere  than  high  values  for  RV  .   Since  the  magnitude 

of  EV^  would  probably  depend  mostly  on  the  pressure  drcp   at 

the   gun   valve,   which  seals  the  breech  chamber,  its  design 

should   te   closely   watched   to   ensure   rapid  ram  times. 

Likewise   the   system   designer   will   have  to   pay  close 

attention  to  the  design  of  the  piping  and  valves  in  a  large 
scale  LPG  feed  system  to  achieve  optimum  performance. 


The  value  cf  the  effective  bulk  modulus  was  determined 
from  the  undamped  natural  freguency  of  the  "water  hammer" 
pressure  oscillations  as  will  be  discussed  in  the  next 
section.  The  bulk  modulus  is  a  fluid  parameter  which 
characterizes  the  spring  effect  of  a  liquid.  The  bulk 
modulus  can  be  substantially  lowered  by  the  elasticity  of 
the  chamber  and  connecting  line  walls  and  the  amcunt  of 
entrapped   gas   present  in  the  propellant  fluid.   As  can  be 
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seen  from  Figure  (11)  ,  the  value  of  the  bulk  modulus  has 
only  a  small  effect  on  ram  time.  Between  bulk  modulus 
values  of  10,000  and  100,000  psi  (the  expected  operating 
region  of  an  operational  system)  the  ram  time  is  almost 
constant,  varying  less  than  a  quarter  of  a  millisecond. 
This  is  an  encouraging  result  since  minimizing  entrapped  gas 
is  a  difficult  and  cpstly  design  constraint. 
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PARAMETEES  AFFECTING  CHAMBER  PRESSURE  OSCILLATIONS 


As  explained  in  Section  III,  after  the  projectile  has 
stepped  the  chamber  pressure  can  he  described  by  a  second 
order  differential  equation  (Egn.  10) .  If  the  system 
described  by  this  differential  equation  is  underdamped,  the 
solution  to  the  equation  is  an  exponentially  damped  sinusoid 
with  a  characteristic  damping  ratio  and  undamped  natural 
freguency. 

From  Egn.  (10),  the  system  damping  ratio,  o   is: 


tf  -5a 


a2  isfpl 


If  this  daaping  ratio  is  increased,  the  peak  chamber 
pressure  as  sell  as  the  dissipation  period  of  the  system 
pressure  oscillations  will  be  decreased.  From  a  designer1 s 
point  of  view,  it  would  be  advantageous  for  both  of  these 
quantities  to  be  as  low  as  possible.  Hence,  it  would  seem 
that  an  increase  in  system  resistance  to  enhance  damping 
would  be  desirable.  However,  as  shown  in  the  previous 
section, an  increase  in  the  system  resistance  also  increases 
the  ram  time-  an  undesirable  occurrence.  Thus,  alterations 
in  system  configuration  which  affect  resistance  coefficient 
values  will  have  to  be  accomplished  carefully  tc  insure 
optimal  system  performance. 

From   Egn.   (1C),   the  undamped  natural  frequency  cf  the 
pressure  cscillations  in  the  LFG  system  is: 

wn  -Vrfr 

T/>L3L3 

The   primary  effect  of  ullage  would  be  to  decrease  the  value 
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of  effective  tulk  modulus  and  consequently  tc  decrease  the 
natural  frequency  of  pressure  oscillations.  Hence,  should 
it  be  found  in  future  tests  that  the  actual  liquid 
prcpellant  is  highly  sensitive  to  pressure  oscillations, 
changes  in  ullage  may  be  effected  to  decrease  system 
ringing.  Fortunately,  as  shown  in  the  previous  section,  the 
effect  of  ullage  variations  on  ram  time  is  minimal. 
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VII.   CONCLGSIONS 


Several  important  conclusions  can  be  drawn  frcm  the 
output  of  the  IPG  ram  feed  computer  nodel.  Thus,  it  is 
realized  that  serious  design  consideration  must  be  given  to 
the  values  cf  fluid  resistances  in  the  feed  system  in  order 
tc  minimize  ram  time  and  peak  pressures,  and  to  optimize  the 
damping  cf  system  pressure  transients.  Fortuantely, 
ccncommitant  model  results  also  show  that  the  effect  of 
entrapped  gas  is  ninimal  for  the  overall  operation  of  a  LPG 
feed  system.  The  model  was  unable  to  account  fcr  any 
cavitation  effects  so  this  is  an  area  that  must  be  a  subject 
cf  further  studies. 

The  fact  that  large  pressure  transients  cccur  during  the 
projectile  ramming  cycle  indicates  the  difficulties  that 
might  arise  with  proposed  designs  for  liquid  propellant  guns 
with  variable  propellant  volumes  and  projectile 
displacements.  The  pressure  transients  experienced  in  these 
systems  will  undcuttedly  be  very  complex  in  nature. 

If  further  improvements  are  desired  in  the  modeling  of 
LPG  feed  systems,  the  next  step  should  be  the  inclusion  of 
thermodynamic  effects.  For  any  real  fluid,  resultant  feed 
system  pressure  and  velocity  changes  will  cause  temperature 
changes  which  may  substantially  affect  the  propellant  fluid 
density.  Consequently,  propellant  ignition  characteristics 
at  high  rates  cf  fire  will  undoubtedly  be  affected. 

The  lumped  parameter  LPG  feed  system  computer  model  has 
been  shown  tc  te  in  favorable  agreement  with  experinental 
data.    It   is   felt   that   this  model  has  been  sufficiently 
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validiated  tc  allow  its  use  in  more  complicated  LSG  feed 
system  designs.  It  is  expected  that  this  model  will  be  used 
by  the  Naval  Ordnance  Station,  Indian  Head,  MarylaDd  to 
assist  in  studying  the  fluid  dynamic  characteristics  of  a  30 
mm  scale  model  1PG  feed  system  becoming  operational  in  June 
1S76.  It  will  be  necessary  to  make  the  obvious  changes  to 
describe  the  different  geometry  of  the  30  mm  scale  uodel. 
The  additional  line  lengths,  and  the  larger  pressure  drops 
due  to  complex  valve  arrangements  in  the  30  mm  system  will 
reguire  different  values  for  model  back  pressure  and 
resistance  coefficients.  The  presence  of  an  accumulator 
near  the  breech  chamber  will  affect  the  value  of  the  system 
camping  ratic.  It  is  possible  that  the  different  geometry 
and  mass  characteristics  of  the  30  mm  model  could  result  in 
the  injector  piston  being  slower  than  the  projectile. 
However  by  incorporating  all  of  these  alterations,  the 
manner  in  which  each  of  the  system  parameters  affect  total 
system  performance  will  be  able  tc  be  predicted  by  the 
versatile  lumped  parameter  computer  model  developed  in  this 
study. 

As  such,  computer  simulation  is  a  useful  design  tool. 
However,  tc  remain  useful,  it  must  be  viewed  in  its  proper 
perspective.  Full  scale  prototype  testing  is  the  only 
conclusive  method  of  demonstrating  ordnance  system 
performance.  unfortunately,  testing  is  both  costly  and  time 
consuming.  Ihus,  computer  simulation,  no  matter  how  simple 
or  complei,  can  be  used  to  designate  critical  testing 
instances  and  to  identify  those  areas  where  design  efforts 
will  be  most  productive.  In  this  way,  the  time  involved  in 
optimizing  system  performance  as  well  as  subseguent 
production  costs  can  be  reduced.  In  explosives  research  and 
ordnance  design,  which  is  still  an  empirical  and  necessarily 
hazardous  science,  computer  simulation  can  be  particularly 
useful  in  specifying  the  correct  approach  for  design  testing 
practices.   In  this   respect,   such   models   as  the   lumped 
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parameter  model  developed  in  this  study  can  fce  considered  as 
critical  signposts  at  some  of  the  many  crossroads  in  systems 
development.  However,  it  must  be  realized  that  as  such  they 
cnly  point  the  way.  The  vehicle  for  arriving  at  an 
operational  system  can  only  be  diligent  research  and 
engineering  based  on  a  measured  progression  of  system 
demonstrations  overseen  by  dedicated  project  managers  aware 
of  the  many  pitfalls  along  the  way. 
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APPENDIX  A 


COMPUTER  PROGRAM  LISTING 


The  fcasis  of  the  state  variable  method  of  systems 
analysis  is  the  interpretation  that  the  energy  state  of  the 
components  of  a  system  completely  describes  the  condition  of 
a  system.  As  defined  in  Ref.  6,  the  state  of  a  system  is  the 
set  of  variables,  the  "state  variables",  which  contain 
sufficient  icformation  about  the  present  condition  of  the 
system  tc  permit  the  determination  of  all  future  time 
history  cf  the  system  -  provided  that  all  future  inputs  to 
the  system  are  known.  Therefore  the  energy  state  cf  those 
elements  which  store  energy  completely  describes  the  system. 
In  the  case  of  the  lumped  parameter  model,  these  elements 
are  the  fluid  inertia  and  capacitance  and  the  inertia  cf  the 
injector  pistcn  and  projectile.  The  energy  state  cf  these 
storage  elements  can  be  described  as  a  function  of  tine  in 
terms  of  the  state  variables  -  pressure,  velocity,  and  the 
system  input  pressure.  Only  those  variables  required  to 
completely  specify  the  state  of  the  system  need  to  be 
included. 

The  first  step  in  arriving  at  the  computer  program  was 
to  take  eguaticns  (1>  through  (9)  and  normalize,  or  solve 
them  for  the  highest  derivative.  The  state  variables  then 
become  the  pressure  or  velocity  associated  with  these  first 
derivatives.  The  computer  program  was  developed  by 
ccmtining  eguations  and  defining  new  variables  as  follows: 
X(1)  =  17 
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X<2)  =  E 

X(3)  =  ?_ 

X<4)  =  ^ 

X(5)  =  1 

1(8)  -  Eg 

X(9)  =  ffi 
Eguation  (1)  becomes: 

-Kfp      At      Ar     At  n 
XDOT(l)  =  — -  X(l)  -  --   X(2)  +  --  Pr  -  --  Pdi 

Mp       Mp      Mp     Mp 

In  the  computer  program,  PR  is  represented  by  X(30)  acd  KFP 

and   PDI   beccme   program   constants,    C  (7)    and   C (2) , 

respectively.   The  value  of  C (2)  and  C(7)  are  input  od  data 

cards.   Dsing  the  fact  that  P  =P  ,  eguation  (2)  becomes: 

1   2 

XD0TC2)  =  p  X(l)  -  -  X(3) 

where   K   is   input   to  the   computer   program  as  C  (3) .  By 

combining  equations  ,(3)  and  (4)  and  using  the  incompressible 

flew  assuopticn  A  V  =A  V  tire  next  eguation  is  derived: 
12   2  4 

XD0K3)  =  -i—  (  X(2)  -  Rv  -i  X(3)  -  X(4)  )    . 
Rho  Lt  A2 

where   BV    tecomes  C  (4) .    Next  eguation   (5)    and   the 

incompressible   flow   assumptions   A  V  =A  V    and   A  V  =A  V 

1224  2536 

yield: 

XDOTW)  =  5l  1  X(3)  -  ^  r  x«) 

A2  l2  a2  l2 

and   as   before,  C  (3) =K.   Eguations  (6)  and  (7)  combine  with 

the  flow  velocity  eguation  A  V  =A  V  to  yield: 

2  5   3  6 

XD°T(5)  =  5^  (  XW  -   RV2  A^  X(5)  -  X(8)  > 

where   El?   becomes  C  (5)  .   Since  V  =V  ,    eguation  (8)  becomes: 
2  6   7 


XD0TC8)  =  --  X(5)  -  r  X(9) 

l_3        L3 


47 


where  K  agair  becomes  C(3) .   Finally,  equation  (9)   becomes: 

-Kfs       A       A 

XDOTO)  ■ X(9)  +  -2  X(8)  -  ?  Pds 

Ms      Ms   .    Ms 

where   KFS   and   EDS  become  C  (8)  and  C  (9)  respectively.   The 
auxiliary  equations: 

XD0T(6)  =  X(l)  ' 
XD0K7)  =  X(9) 

are   used   to  calculate  injector  piston  displacement,  X  (6) , 
and  projectile  displacement,  X (7) . 
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